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Abstract. In this paper, we develop a new censored quantile instrumental variable (CQIV) 
estimator and describe its properties and computation. The CQIV estimator combines Pow- 
ell (1986) censored quantile regression (CQR) to deal semiparametrically with censoring, 
with a control variable approach to incorporate endogenous regressors. The CQIV estima- 
tor is obtained in two stages that are nonadditive in the unobservables. The first stage 
estimates a nonadditive model with infinite dimensional parameters for the control variable, 
such as a quantile or distribution regression model. The second stage estimates a nonaddi- 
tive censored quantile regression model for the response variable of interest, including the 
estimated control variable to deal with endogeneity. For computation, we extend the algo- 
rithm for CQR developed by Chernozhukov and Hong (2002) to incorporate the estimation 
of the control variable. We give generic regularity conditions for asymptotic normality of the 
CQIV estimator and for the validity of resampling methods to approximate its asymptotic 
distribution. We verify these conditions for quantile and distribution regression estimation 
of the control variable. We illustrate the computation and applicability of the CQIV esti- 
mator with numerical examples and an empirical application on estimation of Engel curves 
for alcohol. 
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1. Introduction 

Censoring and endogeneity are common problems in data analysis. For example, income 
survey data are often top-coded and many economic variables such as hours worked, wages 
and expenditure shares are naturally bounded from below by zero. Endogeneity is also an 
ubiquitous phenomenon both in experimental studies due to partial noncompliance (Angrist, 
Imbens, and Rubin, 1996), and in observational studies due to simultaneity (Koopmans 
and Hood, 1953), measurement error (Frish, 1934), sample selection (Heckman, 1979) or 
more generally to the presence of relevant omitted variables. Censoring and endogeneity 
often come together. Thus, for example, we motivate our analysis with the estimation of 
Engel curves for alcohol - the relationship between the share of expenditure on alcohol 
and the household's budget. For this commodity, more than 15% of the households in our 
sample report zero expenditure, and economic theory suggests that total expenditure and 
its composition are jointly determined in the consumption decision of the household. Either 
censoring or endogeneity lead to inconsistency of traditional mean and quantile regression 
estimators by inducing correlation between regressors and error terms. We introduce a 
quantile regression estimator that deals with both problems and name this estimator the 
censored quantile instrumental variable (CQIV) estimator. 

Our procedure deals with censoring semiparametrically through the conditional quantile 
function following Powell (1986). This approach avoids the strong parametric assumptions 
of traditional Tobit estimators. The key ingredient here is the equivariance property of 
quantile functions to monotone transformations such as censoring. Powell's censored quan- 
tile regression estimator, however, has proven to be difficult to compute. We address this 
problem using the computationally attractive algorithm of Chernozhukov and Hong (2002). 
An additional advantage of focusing on the conditional quantile function is that we can cap- 
ture nonadditive heterogeneity in the effects of the regressors across the distribution of the 
response variable by computing CQIV at different quantiles (Koenker, 2005). The traditional 
Tobit framework rules out this heterogeneity by imposing a location shift model. 

We deal with endogeneity using a control variable approach. The basic idea is to add a 
variable to the regression such that, once we condition on this variable, regressors and error 
terms become independent. This so-called control variable is usually unobservable and need 
to be estimated in a first stage. Our main contribution here is to allow for semiparametric 
models with infinite dimensional parameters and nonadditive error terms, such as quantile 
regression and distribution regression, to model and estimate the first stage and back out the 
control variable. This part of the analysis constitutes the main theoretical difficulty because 
the first stage estimators do not live in spaces with nice entropic properties, unlike, for 
example, in Andrews (1994) or Newey (1994). To overcome this problem, we develop a new 
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technique to derive asymptotic theory for two-stage procedures with plugged-in first stage 
estimators that, while not hving in Donsker spaces themselves, can be suitably approximated 
by random functions that live in Donsker spaces. The CQIV estimator is therefore obtained 
in two stages that are nonadditive in the unobservables. The first stage estimates the control 
variable, whereas the second stage estimates a nonadditive censored quantile regression model 
for the response variable of interest, including the estimated control variable to deal with 
endogeneity. 

We analyze the theoretical properties of the CQIV estimator in large samples. Under suit- 
able regularity conditions, CQIV is y^-consistent and has a normal limiting distribution. 
We characterize the expression of the asymptotic variance. Although this expression can be 
estimated using standard methods, we find it more convenient to use resampling methods 
for inference. We focus on weighted bootstrap because it has practical advantages over non- 
parametric bootstrap to deal with discrete regressors with small cell sizes (Ma and Kosorok, 
2005, and Chen and Pouzo, 2009). We give regularity conditions for the consistency of 
weighted bootstrap to approximate the distribution of the CQIV estimator. For our leading 
cases of quantile and distribution regression estimation of the control variable, we provide 
more primitive assumptions that verify the regularity conditions for asymptotic normality 
and weighted bootstrap consistency. The verification of these conditions for quantile and 
distribution regression estimators of the first stage is new to the best of our knowledge. 

The CQIV estimator is simple to compute using standard statistical software. We demon- 
strate its implementation through Monte-Carlo simulations and an empirical application to 
the estimation of Engel curves for alcohol. The results of the Monte- Carlo exercise demon- 
strate that the performance of CQIV is comparable to that of Tobit IV in data generated 
to satisfy the Tobit IV assumptions, and it outperforms Tobit IV under heteroskedasticity. 
The results of the application to Engel curves demonstrate the importance of accounting for 
endogeneity and censoring in real data. Another application of our CQIV estimator to the 
estimation of the price elasticity of expenditure on medical care appears in Kowalski (2009). 

1.1. Literature review. There is an extensive previous literature on the control variable 
approach to deal with endogeneity in models without censoring. Hausman (1978) and 
Wooldridge (2010) discussed parametric triangular linear and nonlinear models. Newey, 
Powell, and Vella (1999) described the use of this approach in nonparametric triangular sys- 
tems of equations for the conditional mean, but limited the analysis to models with additive 
errors both in the first and the second stage. Lee (2007) set forth an estimation strategy 
using a control variable approach for a triangular system of equations for conditional quan- 
tiles with an additive nonparametric first stage. Imbens and Newey (2002, 2009) extended 
the analysis to triangular nonseparable models with nonadditive error terms in both the 
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first and second stage. They focused on identification and nonparametric estimation rates 
for average, quantile and policy effects. Our paper complements Imbens and Newey (2002, 
2009) by providing inference methods and allowing for censoring. Chesher (2003) and Jun 
(2009) considered local identification and semiparametric estimation of uncensored trian- 
gular quantile regression models with a nonseparable control variable. Relative to CQIV, 
these local methods impose less structure in the model at the cost of slower rates of con- 
vergence in estimation. While the previous papers focused on triangular models, Blundell 
and Matzkin (2010) have recently derived conditions for the existence of control variables 
in nonseparable simultaneous equations models. We refer also to Matzkin (2007) for an 
excellent comprehensive review of results on nonparametric identification of triangular and 
simultaneous equations models. 

Our work is also closely related to Ma and Koenker (2006). They considered identifica- 
tion and estimation of quantile effects without censoring using a parametric control variable. 
Their parametric assumptions rule out the use of nonadditive models with infinite dimen- 
sional parameters in the first stage, such as quantile and distribution regression models in 
the first stage. In contrast, our approach is specifically designed to handle the latter, and 
in doing so, it puts the first stage and second stage models on the equally fiexible foot- 
ing. Allowing for a nonadditive infinite dimensional control variable makes the analysis of 
the asymptotic properties of our estimator very delicate and requires developing new proof 
techniques. In particular, we need to deal with control variable estimators depending on 
random functions that do not live in Donsker classes. We address this difficulty approx- 
imating these functions with sufficient degree of accuracy by smoother functions that live 
in Donsker classes. In the case of quantile and distribution regression, we carry out this 
approximation by smoothing the empirical quantile regression and distribution regression 
processes using third order kernels. 

For models with censoring, the literature is more sparse. Smith and Blundell (1986) pio- 
neered the use of the control variable approach to estimate a triangular parametric additive 
model for the conditional mean. More recently, Blundell and Powell (2007) proposed an 
alternative censored quantile instrumental variable estimator that assumes additive errors in 
the first stage. Our estimator allows for a more fiexible nonadditive first stage specification. 



1.2. Plan of the paper. The rest of the paper is organized as follows. In Section |2l we 
present the CQIV model, and develop estimation and inference methods for the parameters of 
interest of this model. In Section [HI we describe the associated computational algorithms and 
present results from a Monte-Carlo simulation exercise. In Section HI we present an empirical 
application of CQIV to Engel curves. In Section [5l we provide conclusions and discuss 
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potential empirical applications of CQIV. The proofs of the main results, and additional 
details on the computational algorithms and numerical examples are given in the Appendix. 

2. Censored Quantile Instrumental Variable Regression 

2.1. The Model. We consider the following triangular system of quantile equations: 

Y = max(F*,C), (2.1) 
Y* = Qy4U\D,W,V), (2.2) 
D = QDiV\W,Z). (2.3) 

In this system, Y* is a continuous latent response variable, the observed variable Y is ob- 
tained by censoring Y* from below at the level determined by the variable C, D is the 
continuous regressor of interest, is a vector of covariates, possibly containing V is a, 
latent unobserved regressor that accounts for the possible endogeneity of D, and Z is a vec- 
tor of "instrumental variables" excluded from fl2.2p FI Further, u Qy*{u \ D, W, V) is the 
conditional quantile function of Y* given {D, W, V); and v Qd{v \ W, Z) is the conditional 
quantile function of the regressor D given {W, Z). Here, [/ is a Skorohod disturbance for Y 
that satisfies the independence assumption 

U ~ U{0,1) I D,W,Z,V,C, 

and is a Skorohod disturbance for D that satisfies 

V ~ U{0, 1) I W, Z, C. 

In the last two equations, we make the assumption that the censoring variable C is inde- 
pendent of the disturbances U and V . This variable can, in principle, be included in W . To 
recover the conditional quantile function of the latent response variable in equation ( 12. 2p . 
it is important to condition on an unobserved regressor V which plays the role of a "con- 
trol variable." Equation (12.31) allows us to recover this unobserved regressor as a residual 
that explains movements in the variable -D, conditional on the set of instruments and other 
covariates. 

In the Engel curve application, Y is the expenditure share in alcohol, bounded from below 
at C = 0, Z) is total expenditure on nondurables and services, W are household demographic 
characteristics, and Z is labor income measured by the earnings of the head of the household. 
Total expenditure is likely to be jointly determined with the budget composition in the 
household's allocation of income across consumption goods and leisure. Thus, households 

"'^We focus on left censored response variables without loss of generality. If Y is right censored dX C, Y ^ 
min(y*,C), the analysis of the paper applies without change to F = —Y, Y* = —Y*, C = —C, and 
Qy* = —Qy* , because Y = max(F*, C). 
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with a high preference to consume "non-essential" goods such as alcohol tend to expend a 
higher proportion of their incomes and therefore to have a higher expenditure. The control 
variable V in this case is the marginal propensity to consume, measured by the household 
ranking in the conditional distribution of expenditure given labor income and household 
characteristics. This propensity captures unobserved preference variables that affect both the 
level and composition of the budget. Under the conditions for a two stage budgeting decision 
process (Gorman, 1959), where the household first divides income between consumption 
and leisure/labor and then decide the consumption allocation, some sources of income can 
provide plausible exogenous variation with respect to the budget shares. For example, if 
preferences are weakly separable in consumption and leisure/labor, the consumption budget 
shares do not depend on labor income given the consumption expenditure (see, e.g., Deaton 
and Muellbauer, 1980). This justifies the use of labor income as an exclusion restriction. 

An example of a structural model that has the triangular representation f l2.2p - fl2.3p is the 
following system of equations: 

Y* = gy{D,W,e), (2.4) 
D = gniW,Z,V), (2.5) 

where gy and g^ are increasing in their third arguments, and e ~ f/(0, 1) and V ^ U{0, 1) 
independent of (W, Z, C). By the Skorohod representation for e, e = Qe{U \ V) = ge{V, U), 
where U ~ U{0,1) independent of {D,W, Z,V,C). The corresponding conditional quantile 
functions have the form of (12. 2p and (12.30 with 

Qy,{u\D,W,V) = gy{D,W,g,{V,u)), 
Qd{v\W,Z) = gD{W,Z,v). 

In the Engel curve application, we can interpret V as the marginal propensity to consume out 
of labor income and U as the unobserved household preference to spend on alcohol relative 
to households with the same characteristics W and marginal propensity to consume V. 

In the system of equations (I2.ip - (l2.3p . the observed response variable has the quantile 
representation 

Y = Qy{U\ D,W,V,C) = max{Qy*{U \ D,W,V),C), (2.6) 

by the equivariance property of the quantiles to monotone transformations. For example, 
the quantile function for the observed response in the system of equations (I2.4p - (l2.5p has 
the form: 

Qy{u I D, W, V, C) = max{^y (D, W, g,{V, u)), C}. 
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Whether the response of interest is the latent or observed variable depends on the source of 
censoring (e.g., Wooldridge, 2010). When censoring is due to data limitations such as top- 
coding, we are often interested in the conditional quantile function of the latent response 
variable Qy and marginal effects derived from this function. For example, in the system 
( I2.4p -( E3]1 the marginal effect of the endogenous regressor D evaluated at {D, W, V, U) = 
{d, w, V, u) is 

ddQY*{u I d,w,v) = ddgY{d,w,ge{v,u)), 

which corresponds to the ceteris paribus effect of a marginal change of D on the latent 
response Y* for individuals with {D,W,e) = {d,w, g^{v,u)). When the censoring is due 
to economic or behavioral reasons such are corner solutions, we are often interested in the 
conditional quantile function of the observed response variable Qy and marginal effects 
derived from this function. For example, in the system (I2.4p - fl2.5p the marginal effect of the 
endogenous regressor D evaluated at {D, W, V, U, C) = {d, w, v, u, c) is 

ddQviu I d,w,v,c) = l{gY{d,w,g,{v,u)) > c}ddgYid,w, g,{v,u)), 

which corresponds to the ceteris paribus effect of a marginal change of D on the observed 
response Y for individuals with {D, W, e, C) = {d, w, ge{v, u),c). Since either of the marginal 
effects might depend on individual characteristics, average marginal effects or marginal effects 
evaluated at interesting values are often reported. 



2.2. Generic Estimation. To make estimation both practical and realistic, we impose a 
flexible semiparametric restriction on the functional form of the conditional quantile function 
in (12. 2p . In particular, we assume that 

Qy,{u\D,W,V)=X'Po{u), X = x{D,W,V), (2.7) 

where x{D,W,V) is a vector of transformations of the initial regressors {D,W,V). The 
transformations could be, for example, polynomial, trigonometric, B-spline or other basis 
functions that have good approximating properties for economic problems. An important 
property of this functional form is linearity in parameters, which is very convenient for 
computation. The resulting conditional quantile function of the censored random variable 

F = max(F*,C), 

is given by 

Qy{u I D, W, V, C) = max(X'/3o(M), C). (2.8) 

This is the standard functional form for the censored quantile regression (CQR) first derived 
by Powell (1984) in the exogenous case. 
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Given a random sample {Yi, Di,Wi, Zi,Ci}^^^, we form the estimator for the parameter 

/3o{u) as 

1 " ^ 

/3(M) = arg min -TUS'ff > ^)pu{Yi - X'^P), (2.9) 

;3eRdim(X) n ^-^ 

where Pu{z) = {u — l{z < 0))2; is the asymmetric absolute loss function of Koenker and 
Bassett (1978), Xj = x{Di,Wi,Vi), Si = s(Xj,Cj), s{X,C) is a vector of transformations 
of (X, C), and Vi is an estimator of Vi. This estimator adapts the algorithm for the CQR 
estimator developed in Chernozhukov and Hong (2002) to deal with endogeneity. We call 
the multiplier 1(5*17 > ^) the selector, as its purpose is to predict the subset of individuals 
for which the probability of censoring is sufficiently low to permit using a linear - in place 
of a censored linear - functional form for the conditional quantile. We formally state the 
conditions on the selector in the next subsection. The estimator in (12. 9p may be seen as a 
computationally attractive approximation to Powell estimator applied to our case: 

1 " 

/3p{u) = arg min - V^p^fFi - max(X-/3, Ci)]. 

ag|5dim(X) n ^ 

The CQIV estimator will be computed using an iterative procedure where each step will 
take the form specified in equation (12. 9p . We start selecting the set of "quantile- uncensored" 
observations for which the conditional quantile function is above the censoring point. We 
implement this step by estimating the conditional probabilities of censoring using a fiexible 
binary choice model. Quantile-uncensored observations have probability of censoring lower 
than the quantile index u. We estimate the linear part of the conditional quantile function, 
Xj'/3o(m), on the sample of quantile-uncensored observations by standard quantile regression. 
Then, we update the set of quantile-uncensored observations by selecting those observations 
with conditional quantile estimates that are above their censoring points and iterate. We 
provide more practical implementation details in the next section. 

The control variable V can be estimated in several ways. Note that if Qd{v \ W, Z) is 
invertible in v, the control variable has several equivalent representations: 

V = MD, W, Z) = Fd{D \W,Z) = Q-^\D \W,Z)= [ 1{Qd{v \W,Z)< D]dv. (2.10) 

Jo 

For any estimator of Fd{D \ W, Z) or Qd{V \W,Z), denoted by Fd(£' | W, Z) or Qd{V \ 
W, Z), based on any parametric or semi-parametric functional form, the resulting estimator 
for the control variable is 

V = d{D, W, Z) = Fd{D \W,Z) or V = d{D, W, Z) = [ liQoiv \ W, Z) < D}dv. 

Jo 

Here we consider several examples: in the classical additive location model, we have that 
Qniv \W, Z) = R'ttq + Qv{v), where Qy is a quantile function, and R = r(W, Z) is a vector 
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collecting transformations of W and Z. The control variable is 

V = Qy\D-R'7ro), 

which can be estimated by the empirical CDF of the least squares residuals. Chernozhukov, 
Fernandez- Val and Melly (2009) developed asymptotic theory for this estimator. If D \ 
W, Z A^(i?'7ro, 0"^), the control variable has the common parametric form V = $~^([Z) — 
R'TTol/a), where denotes the quantile function of the standard normal distribution. This 
control variable can be estimated by plugging in estimates of the regression coefficients and 
residual variance. 

In a non-additive quantile regression model, we have that Qd{v \ W, Z) = R'ttq^v), and 

V = Qd\D \W,Z)= [ l{R'no{v) < D}dv. 

Jo 

The estimator takes the form 

[ l{R'7r{v) < D}dv, (2.11) 
Jo 

where n{v) is the Koenker and Bassett (1978) quantile regression estimator and the integral 
can be approximated numerically using a finite grid of quantiles. The use of the integral to 
obtain a generalized inverse is convenient to avoid monotonicity problems in v i— )■ R'ti{v) due 
to misspecification or sampling error. Chernozhukov, Fernandez- Val, and Galichon (2010) 
developed asymptotic theory for this estimator. 

We can also estimate ■j?o using distribution regression. In this case we consider a semi- 
parametric model for the conditional distribution of D to construct a control variable 

V = Fr>{D I 1^,Z) =A(i?Vo(D)), 

where A is a probit or logit link function. The estimator takes the form 

V ^ K{R'tv{D)), 

where 7r((i) is the maximum likelihood estimator of 7ro((i) at each d (see, e.g., Foresi and Per- 
acchi, 1995, and Chernozhukov, Fernandez- Val and Melly, 2009). Chernozhukov, Fernandez- 
Val and Melly (2009) developed asymptotic theory for this estimator. 

2.3. Regularity Conditions for Estimation. In what follows, we shall use the following 

notation. We let the random vector A = {Y, D,W, Z,C, X,V) live on some probability 
space {Qo, J-'o, P). Thus, the probability measure P determines the law of A or any of its 
elements. We also let Ai,...,An, i.i.d. copies of A, live on the complete probability space 
(Q, J-", P), which contains the infinite product of (Qq, J^o, P)- Moreover, this probability space 
can be suitably enriched to carry also the random weights that will appear in the weighted 
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bootstrap. The distinction between the two laws P and P is helpful to simplify the notation 
in the proofs and in the analysis. Calligraphic letter such as 3^ and X denote the support 
of Y and X; and yX denotes the joint support of {Y,X). Unless explicitly mentioned, all 
functions appearing in the statements are assumed to be measurable. 

We now state formally the assumptions. The first assumption is our model. 

Assumption 1 (Model). We have {¥{, Di,Wi, Zi,Ci}^^i, a sample of size n of independent 
and identically distributed observations from the random vector (Y, D, W, Z, C) which obeys 
the model assumptions stated in equations ^2. l\j - Ii2.10\) . i.e. 

Qy{u I D, W, Z, V, C) =Qy{u\ X, C) = max(X'/3o(M), C), X = x{D, W, V), 
V = MD, W, Z) = Fd{D \ W,Z)^ f/(0, 1) I W, Z. 

The second assumption imposes compactness and smoothness conditions. Compactness 
can be relaxed at the cost of more complicated and cumbersome proofs, while the smoothness 
conditions are fairly tight. 

Assumption 2 (Compactness and smoothness), (a) The set yVWZCX is compact, (b) 
The endogenous regressor D has a continuous conditional density //)(■ | w, z) that is bounded 
above by a constant uniformly in {w, z) G WZ. (c) The random variable Y has a conditional 
density fyiy \ x, c) on (c, oo) that is uniformly continuous iny E (c, oo) uniformly in {x, c) G 
XC, and bounded above by a constant uniformly in {x,c) G XC. (d) The derivative vector 
dvx{d,w,v) exists and its components are uniformly continuous in v E [0,1] uniformly in 
{d, w) G VyV, and are bounded in absolute value by a constant, uniformly in {d, w, v) G 

vwv. 

The following assumption is a high-level condition on the function-valued estimator of 
the control variable. We assume that it has an asymptotic functional linear representation. 
Moreover, this functional estimator, while not necessarily living in a Donsker class, can be 
approximated by a random function that does live in a Donsker class. We will fully verify this 
condition for the case of quantile regression and distribution regression under more primitive 
conditions. 

Assumption 3 (Estimator of the control variable). We have an estimator of the control 
variable of the form V = ^{D, W, Z), such that uniformly over {d, w, z) G VWZ, (a) 

1 " 

y/n{'d{d, w, z) - i3o{d, w, z)) = i{Ai, d, w, z) + op(l), Ep[£(A, d, w, z)] = 0, 

where Ep[i{A,D,W,Zf]<oo and ||^ ELi^l^^- OIU = Op(l), and (b) 

||^-^||oo = op(l/v^), for deT, 
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where the entropy of the function class T is not too high, namely 

logiV(e, T, II • lU) < l/(elog2(l/e)), for allO<e<l. 

The following assumptions are on the selector. The first part is a high-level condition on 
the estimator of the selector. The second part is a smoothness condition on the index that 
defines the selector. We shall verify that the CQIV estimator can act as a legitimate selector 
itself. Although the statement is involved, this condition can be easily satisfied as explained 
below. 

Assumption 4 (Selector), (a) The selection rule has the form 

l[s{x{D,W,V),Cyj>c;], 

for some <^ > 0, where 7 — J-p 70 and, for some e' > 0, 

l[S'7o > ^/2] < l[X'/3o{u) >C + e']< 1|X'/3o(m) > C] P-a.e., 

where S = s{X,V) and 1[X'/3o(m) > C] = 1[P{Y = C \ Z,W,V) < u]. (h) The set S 
is compact, (c) The density of the random variable s{x{D,W,'d{D,W, Z)),C)''~^ exists and 
is hounded above by a constant, uniformly in E T and in d G T, where T is an open 
neighborhood of '-/o and T is defined in Assumption\^ (d) The components of the derivative 
vector djjs{x{d, w,v), a) are uniformly continuous at each v G [0,1], uniformly in {d,w,c) G 
VyVC, and are bounded in absolute value by a constant, uniformly in {d,w,v,c) G VWVC. 

The next assumption is a sufficient condition to guarantee local identification of the pa- 
rameter of interest as well as A/n-consistency and asymptotic normality of the estimator. 

Assumption 5 (Identification and non- degeneracy) . (a) The matrix 

J{u) := Ep[/y(X%(M) I X,C)XX' l(5'7o > ^)] 

is of full rank, (b) The matrix 

A{u):=VaTp[f{A)+g{A)], 
is finite and is of full rank, where 

fiA) := {!{¥ < X'/3o(n)) - n}Xl(5'7o > 0, 
and, for X = d^x^D, W,v)\v=v! 

g{A) := Ep[fy{X'f3o{u) \ X, C)XX%(m)1(5'7o > ^^^a, D,W, Z)]l^^. 

Assumption 4(a) requires the selector to find a subset of the quantile-censored observations, 
whereas Assumption 5 requires the selector to find a nonempty subset. Given /3°(m), an initial 
consistent estimator of /3o{u), we can form the selector as l[s{x{D, W, V), C)'7 > where 
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s{x{D,W,V),C) = [x{D,W,Vy,C]', 7 = [^^u)',-!]', and ^ is a small fixed cut-off that 
ensures that the selector is asymptotically conservative but nontrivial. To find we use 

a selector based on a flexible model for the probability of censoring. This model does not 
need to be correctly specified under a mild separating hyperplane condition for the quantile- 
uncensored observations (Chernozhukov and Hong, 2002). Alternatively, we can estimate a 
fully nonparametric model for the censoring probabilities. We do not pursue this approach 
to preserve the computational appeal of the CQIV estimator. 

2.4. Main Estimation Results. The following result states that the CQIV estimator is 
consistent, converges to the true parameter at a ^/n rate, and is normally distributed in large 
samples. 

Theorem 1 (Asymptotic distribution of CQIV). Under the stated assumptions 

MMu) - f3o{u)) N{0, J-\u)A{u)J-\u)). 

We can estimate the variance-covariance matrix using standard methods and carry out 
analytical inference based on the normal distribution. Estimators for the components of the 
variance can be formed, e.g., following Powell (1991) and Koenker (2005). However, this 
is not very convenient for practice due to the complicated form of these components and 
the need to estimate conditional densities. Instead, we suggest using weighted bootstrap 
(Chamberlain and Imbens, 2003, Ma and Kosorok, 2005, Chen and Pouzo, 2009) and prove 
its validity in what follows. 

We focus on weighted bootstrap because it has practical advantages over nonparametric 
bootstrap to deal with discrete regressors with small cell sizes and the proof of its consistency 
is not overly complex, following the strategy set forth by Ma and Kosorok (2005). Moreover, 
a particular version of the weighted bootstrap, with exponentials acting as weights, has a 
nice Bayesian interpretation (Chamberlain and Imbens, 2003). 

To describe the weighted bootstrap procedure in our setting, we first introduce the "weights" . 

Assumption 6 (Bootstrap weights). The weights (ei, e„) are i.i.d. draws from a random 
variable e > 0, with Ep[e] = 1 and Varp[e] = 1, living on the probability space {Q,J-',¥) and 
are independent of the data {Yi, Di, Wi, Zi, Cijf^i for all n. 

Remcirk 1 (Bootstrap weights). The chief and recommended example of bootstrap weights 
is given by e set to be the standard exponential random variable. Note that for other 
positive random variables with Ep[e] = 1 but Varp[e] > 1, we can take the transformation 
e — l + {e — 1)/Varp[e]^/^, which satisfies e > 0, Ep[e] = 1, and Varp[e] = 1. 

The weights act as sampling weights in the bootstrap procedure. In each repetition, we 
draw a new set of weights (ei, . . . , e„) and recompute the CQIV estimator in the weighted 
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sample. We refer to the next section for practical details, and here we define the quantities 
needed to verify the vahdity of this bootstrap scheme. Specifically, let denote the esti- 
mator of the control variable for observation i in the weighted sample, such as the quantile 
regression or distribution regression based estimators described in the next section. The 
CQIV estimator in the weighted sample solves 

1 " ^ 

= arg imn -J2e^l{lSt><;)Pu{Y,-(3'X^), (2.12) 

where Xf = x(Di,Wi,Vf). = s(Xf,Cj), and 7 is a consistent estimator of the selector. 
Note that we do not need to recompute 7 in the weighted samples, which is convenient for 
computation. 

We make the following assumptions about the estimator of the control variable in the 
weighted sample. 

Assumption 7 (Weighted estimator of control variable). Let (ei, . . . , e„) be a sequence of 
weights that satisfies Assumption 6. We have an estimator of the control variable of the form 
ye _ ^e^]j^ yy^ z), such that uniformly over VVVZ, 

1 

\/n09'^id, w, z) — doid, w, z)) = —;= > eMAi, d, w, z) + op(l), Ep[£(A, d, w, z)] = 0, 

where'E.p[£{A,D,W,Zf]< 00 and ||^ ELi e^A, OIU = Op(l), and 

||^'-^1loo = op(l/V^), for ^'eT, 
where the entropy of the function class T is not too high, namely 

\ogN{e, T, II • \U < l/(elog2(l/e)), for a// < e < 1. 

Basically this is the same condition as Assumption 3 in the unweighted sample, and 
therefore both can be verified using analogous arguments. Note also that the condition is 
stated under the probability measure P, i.e. unconditionally on the data, which actually 
simplifies verification. We give primitive conditions that verify this assumption for quantile 
and distribution regression estimation of the control variable in the next section. 

The following result shows the consistency of weighted bootstrap to approximate the 
asymptotic distribution of the CQIV estimator. 

Theorem 2 (Weighted-bootstrap validity for CQIV). Under the stated assumptions, condi- 
tionally on the data 

V^0%u) - P{u)) N{0, J-\u)A{u)J-\u)), 

in probability under P. 
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Note that the statement above formally means that the distance between the law of 
^/n{/3^{u)—/3{u)) conditional on the data and the law of the normal vector A^(0, J~^{u)A{u)J^ 
as measured by any metric that metrizes weak convergence, conveges in probability to zero. 
More specifically, 

dBL{C[V^{P^{u) - ^{u))\data], C[N{0, J~-\u)A{u)J-\u))]} 0, 

where dsi denotes the bounded Lipshitz metric. 

In practice, we approximate numerically the distribution of ^/n{P'^{u) — (3{u)) conditional 
on the data by simulation. For = 1, . . . , 5, we compute /3^{u) solving the problem fl2.12p 
with the data fixed and a set of weights {eu, e„f,) randomly drawn for a distribution that 
satisfies Assumption 6. By Theorem 2, we can use the empirical distribution of A/n(/3^(u) — 
(3{u)) to make asymptotically valid inference on (3q{u) for large B. 

2.5. Quantile and distribution regression estimation of the control variable. One 

of the main contributions of this paper is to allow for quantile and distribution regression 
estimation of the control variable. The difficulties here are multifold, since the control vari- 
able depends on the infinite dimensional function vro(-), and more importantly the estimated 
version of this function, 7r(-), does not seem to lie in any class with good entropic properties. 
We overcome these difficulties by demonstrating that the estimated function can be approx- 
imated with sufficient degree of accuracy by a random function that lies in a class with good 
entropic properties. To carry out this approximation, we smooth the empirical quantile re- 
gression and distribution regression processes by third order kernels, after suitably extending 
the processes to deal with boundary issues. Such kernels can be obtained by reproducing ker- 
nel Hilbert space methods or via twicing kernel methods (Berlinet, 1993, and Newey, Hsieh, 
and Robins, 2004). In the case of quantile regression, we also use results of the asymptotic 
theory for rearrangement-related operators developed by Chernozhukov, Fernandez- Val and 
Galichon (2010). Moreover, all the previous arguments carry over weighted samples, which 
is relevant for the bootstrap. 

2.5.1. Quantile regression. We impose the following condition: 

Assumption 8 (QR control variable), (a) The conditional quantile function of D given 
{W, Z) follows the quantile regression model, i.e., 

Qd{- \W,Z) = Qd{- I R) = R'M-), R = r{W, Z), 

where the coefficients v i-> 7ro(f ) are three times continuously differentiable with uniformly 
bounded derivatives, and TZ is compact; (b) The conditional density //)(■ | R) is uniformly 
bounded by a constant P-a.e., and is continuous at R'ttq{v) uniformly in v E (0,1) P-a.e. 
(c) The Gram matrix E[RR'] is finite and has full rank. 
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For p^{z) := {y - l{z < 0))z, let 

1 " 

n^{v) e arg min - eip^{Di - R[it), 

«=i 

where either Cj = 1 for the unweighted sample, to obtain the estimates; or Cj is drawn from 
a positive random variable with unit mean and variance for the weighted sample, to obtain 
bootstrap estimates. Then set 

§^{d,r)= \ l{r'7ro{v) <d}dv; d'{d,r) = [ l{r'n^{v) < d}dv. 
i(o,i) i(o,i) 

The following result verifies that our main high-level conditions for the control variable 
estimator in Assumptions 3 and 7 hold under Assumption 8. The verification is done simul- 
taneously for weighted and unweighted samples by including weights that can be equal to 
the trivial unit weights, as mentioned above. 

Theorem 3 (Validity of Assumptions 3 & 7 for QR). Suppose that Assumption\^holds. (1) 
We have that 

1 " 

^{r{d,r)-d^{d,r)) = ed{Ai,d,r) + o^{l) ^ ^\d,r) mr{Vn), 

i{A,d,r) := Md \ r)r'Ep [fDiR'MMd,r)) \ R)RRr' x 
x[l{D< R'TioiMd, r))} - Md, r)]R, 
Ep[i{A,d,r)] = 0, Ep[i{A,D,R)'^] < oo, 

where A'^{d,r) is a Gaussian process with continuous paths and covariance function given 
by Ep[i{A,d,r)i{A,d,ry]. (2) Moreover, there exists : T>TZ H- [0,1] that obeys the same 
first order representation, is close to -d^ in the sense that '^'d^ — ^^\\oo = Of{l/y/n), and, with 
probability approaching one, belongs to a bounded function class T such that 

logAr(e,T,||-|U)<e-i/2, < e < 1. 

Thus, Assumption 3 holds for the case Cj = 1, and Assumption 7 holds for the case of Cj 
being drawn from a positive random variable with unit mean and variance as in Assumption 
6. Thus, the results of Theorem 1 and 2 apply for the QR estimator of the control variable. 

2.5.2. Distribution regression. We impose the following condition: 

Assumption 9 (DR control variable), (a) The conditional distribution function of D given 
{W, Z) follows the distribution regression model, i.e., 

Fn{- \W,Z) = Fn{- \ R) = A{R'rro{-)), R = r{W, Z), 
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where A is either the probit or logit link function, the coefficients d i-> 7ro(d) are three times 
continuously differentiable with uniformly bounded derivatives; (b) V and TZ are compact; 
(c) The Gram matrix ERE! has full rank. 



Let 

1 

e arg min - Vei{l(A < d) \ogK{R'-K) + 1(A > d) log[l - A(i?-7r)]}, 

^glRdim(H) fl ^-^ 
1=1 

where either ti — l for the unweighted sample, to obtain the estimates; or Cj is drawn from 
a positive random variable with unit mean and variance for the weighted sample, to obtain 
bootstrap estimates. Then set 

Md,r) = A(/7ro(d)); d\d,r) = A(r'r(d)). 

The following result verifies that our main high-level conditions for the control variable 
estimator in Assumptions 3 and 7 hold under Assumption 9. The verification is done simul- 
taneously for weighted and unweighted samples by including weights that can be equal to 
the trivial unit weights. 



Theorem 4 (Validity of Assumptions 3 & 7 for DR) . Suppose that Assumption 9 holds. (1) 
We have that 



^/E{d'{d,r)-Md,r)) = 
£{A,d,r) := 

Ep[e(A,d,r)] = 



n 

J2 eAA, d, r) + op(l) w A\d, r) m r{Vn), 
dA{R'no{d)f 



aA(r'7ro((i))r'Ei 



H -1 



X 



A(R'7ro(d))[l - A{R'7ro{d))] 
l{D<d}-HR'Md)) 



RR' 



X 



A(R'7ro{d))[l-A(R'7ro(d))] 
0,Ep[£(A,D,R)^] < oo. 



where A'^{d,r) is a Gaussian process with continuous paths and covariance function given 
by Ep[i{A,d,r)i{A,d,fy], and dA is the derivative of A. (2) Moreover, there exists : 
VTZ I— >■ [0, 1] that obeys the same first order representation, is close to ■Q^ in the sense that 
11''^^ "^^"^ I loo = op(l/\/n) and, with probability approaching one, belongs to a bounded function 
class T such that 

logA^(e,T,||.||oo)<e-V2, < e < 1. 

Thus, Assumption 3 holds for the case ei — 1, and Assumption 7 holds for the case of Cj 
being drawn from a positive random variable with unit mean and variance as in Assumption 
6. Thus, the results of Theorem 1 and 2 apply for the DR estimator of the control variable. 
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3. Computation and Numerical Examples 

This section describes the numerical algorithms to compute the CQIV estimator and 
weighted bootstrap confidence intervals, and shows the results of a Monte Carlo numerical 
example. 

3.1. CQIV Algorithm. The algorithm to obtain CQIV estimates is similar to Chernozhukov 
and Hong (2002). We add an initial step to estimate the control variable V. We name this 
step as Step to facilitate comparison with the Chernozhukov and Hong (2002) 3-Step CQR 
algorithm. 

Algorithm 1 (CQIV). For each desired quantile u, perform the following steps: 

0. Obtain an estimate of the control variable for each individual, Vi, and construct Xi = 

1. Select a subset of quantile-uncensored observations, Jq, whose conditional quantile 
function is likely to be above the censoring point, namely select a subset of {i : 
X-Po{u) > C}. To find these observations, we note that X'Pq{u) > C is equivalent 
to P{Y > C \ X,C) > 1 — u. Hence we predict the quantile-uncensored observations 
using a flexible binary choice model: 

P{Y>C\ X, C) = A{S[5o), = s(X„ C), 

where A is a known link function, typically a probit or a logit. In estimation, we 
replace Si by Si = s(Vj, Cj). Then, we select the sample Jq according to the following 
criterion: 

Jo = {i:A{S'J)>l-u + ko}. 

2. Estimate a standard quantile regression on the subsample defined by Jq: 

^V) = arg min Vp„(V.-V:/3). (3.1) 

le Jo 

Next, using the predicted values, select another subset of quantile-uncensored obser- 
vations, Ji, from the full sample according to the following criterion: 

Ji = {i: X[p\u) >C, + ^i}. (3.2) 

3. Estimate a standard quantile regression on the subsample defined by Ji. Formally, 
replace Jq by Ji in \3.1\) . The new estimates, I3^{u), are the 3-Step CQIV coefficient 
estimates. 

4. (Optional) With the results from the previous step, select a new sample J2 replacing 
/3^{u) by /3^{u) in ^3.2\) . Iterate this and the previous step a bounded number of times. 
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Remark 2 (Step 0). A simple additive strategy is to estimate the control variable using 
the empirical CDF of the residuals from the first stage OLS regression of D on W and Z. 
More flexible non-additive strategies based on quantile regression or distribution regression 
are described in the previous section. 

Remark 3 (Step 1). To predict the quantile-uncensored observations, a probit, logit, or any 
other model that fits the data well can be used. Note that the model does not need to be 
correctly specified; it suffices that it selects a nontrivial subset of observations with X-/3o(m) > 
Ci. To choose the value of fco, it is advisable that a constant fraction of observations satisfying 
A(S'j'5) > 1 — u are excluded from Jq for each quantile. To do so, set fco as the qoth quantile 
of A(S'-5) conditional on A(S'-(5) > 1 — u, where go is a percentage (10% worked well in our 
simulation). The empirical value of ko and the percentage of observations retained in Jq can 
be computed as simple robustness diagnostic tests at each quantile. 

Remark 4 (Step 2). To choose the cut-off ^i, it is advisable that a constant fraction of 
observations satisfying Xl/3^{u) > Ci are excluded from Ji for each quantile. To do so, 
set ^1 to be the gith quantile of XlP^{u) — Ci conditional on Xj'/3°(u) > Ci, where qi is a 
percentage less than go (3% worked well in our simulation). In practice, it is desirable that 
Jo C Ji. If this is not the case, we recommend altering go, gi, or the specification of the 
regression models. At each quantile, the empirical value of <^i, the percentage of observations 
from the full sample retained in Ji, the percentage of observations from Jq retained in Ji, 
and the number of observations in Ji but not in Jo can be computed as simple robustness 
diagnostic tests. The estimator I3^{u) is consistent but will be inefficient relative to the 
estimator obtained in the subsequent step. 

Remark 5 (Steps 1 and 2). In the notation of Assumption HI the selector of Step 1 can be 
expressed as l{S'i^ > <Jo), where S'-7 = Si6 — A^^{l—u) and = {1 — u + k^) — A'^ {1 —u) . 
The selector of Step 2 can also be expressed as l(5'-7 > ^i), where Si = {Xl,Ci)' and 
7 = (^°(«)',-l)'. 

Remark 6 (Steps 2, 3 and 4). Beginning with Step 2, each successive iteration of the 
algorithm should yield estimates that come closer to minimizing the Powell objective func- 
tion. As a simple robustness diagnostic test, we recommend computing the Powell objective 
function using the full sample and the estimated coefficients after each iteration, starting 
with Step 2. This diagnostic test is computationally straightforward because computing the 
objective function for a given set of values is much simpler than maximizing it. In practice, 
this test can be used to determine when to stop the CQIV algorithm for each quantile. If 
the Powell objective function increases from Step s to Step s + 1 for s > 2, estimates from 
Step s can be retained as the coefficient estimates. 
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3.2. Weighted Bootstrap Algorithm. We recommend obtaining confidence intervals through 
a weighted bootstrap procedure, though analytical formulas can also be used. If the esti- 
mation runs quickly on the desired sample, it is straightforward to rerun the entire CQIV 
algorithm B times weighting all the steps by the bootstrap weights. To speed up the com- 
putation, we propose a procedure that uses a one-step CQIV estimator in each bootstrap 
repetition. 

Algorithm 2 (Weighted bootstrap CQIV). For h — 1, . . . ,B , repeat the following steps: 

1. Draw a set of weights (eit, . . . , enb) i-i-d from a random variable e that satisfies As- 
sumption 6. For example, we can draw the weights from a standard exponential 
distribution. 

2. Reestimate the control variable in the weighted sample, = i^l{Di,Wi, Zi), and 
construct X^^ = x{Di, Wi, V^l). 

3. Estimate the weighted quantile regression: 



where Jn, = {i : /3(-u)'Xf5 > Ci + <^i}, and fi{u) is a consistent estimator of Po{u), 
e.g., the 3-stage CQIV estimator {3^ {u). 

Remeirk 7 (Step 2). The estimate of the control function -dl can be obtained by weighted 
least squares, weighted quantile regression, or weighted distribution regression. 

Remark 8 (Step 3). A computationally less expensive alternative is to set Jif, = Ji in all the 
repetitions, where Ji is the subset of selected observations in Step 2 of the CQIV algorithm. 

We can construct an asymptotic (1 — Q;)-confidence interval for a function of the parameter 
vector g{/3o{u)) as [ga/2,9i-a/2\, where ga is the sample a-quantile of [g0l{u)), g0Bi'^))]. 
For example, the 0.025 and 0.975 quantiles of {/3l f.{u), . . . , P% f.{u)) form a 95% asymptotic 
confidence interval for the kth coefficient /5o,fc(i*)- 

3.3. Monte-Carlo illustration. The goal of the following numerical example is to com- 
pare the performance of CQIV relative to tobit IV and other quantile regression estimators 
in finite samples. We generate data according to a normal design that satisfies the tobit 
parametric assumptions and a design with heteroskedasticity in the first stage equation for 
the endogenous regressor D that does not satisfy the tobit parametric assumptions. To 
facilitate the comparison, in both designs we consider a location model for the response vari- 
able Y*, where the coefficients of the conditional expectation function and the conditional 
quantile function are equal (other than the intercept), so that tobit and CQIV estimate the 
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same parameters. A comparison of the dispersion of the tobit estimates to the dispersion 
of the CQIV estimates at each quantile in the first design serves to quantify the relative 
efficiency of CQIV in a case where tobit IV can be expected to perform as well as possible. 
The appendix provides a more detailed description of the designs. 

We consider two tobit estimators for comparison. Tobit-iv is the full information maxi- 
mum likelihood estimator developed by Newey (1987), which is implemented in Stata with 
the command ivtobit. Tobit-cmle is the conditional maximum likelihood tobit estimator 
developed by Smith and Blundell (1986), which uses least squares residuals as a control vari- 
able. For additional comparisons, we present results from the censored quantile regression 
(cqr) estimator of Chernozhukov and Hong (2002), which does not address endogeneity; the 
quantile instrumental variables estimator (qiv-ols) of Lee (2007) with parametric first and 
second stage, which does not account for censoring; and the quantile regression (qr) estima- 
tor of Koenker and Bassett (1978), which does not account for endogeneity nor censoring. 
For CQIV we consider three different methods to estimate the control variable: cqiv-ols, 
which uses least squares; cqiv-qr, which uses quantile regression; and cqiv-dr, which uses 
probit distribution regression. The appendix also provides technical details for all CQIV 
estimators, as well as diagnostic test results for the cqiv-ols estimator. 

We focus on the coefficient on the endogenous regressor D. We report mean bias and root 
mean square error (rmse) for all the estimators at the {.05, .10, .95} quantiles. For the 
homoskedastic design, the bias results are reported in the upper panel of Figure [T] and the 
rmse results are reported in the lower panel. In this figure, we see that tobit-cmle represents 
a substantial improvement over tobit-iv in terms of mean bias and rmse. Even though 
tobit-iv is theoretically efficient in this design, the CQIV estimators out-perform tobit-iv, 
and compare well to tobit-cmle. The figure also demonstrates that the CQIV estimators 
out-perform the other quantile estimators at all estimated quantiles. All of our qualitative 
findings hold when we consider unreported alternative measures of bias and dispersion such 
as median bias, interquartile range, and standard deviation. 

The similar performance of tobit-cmle and cqiv can be explained by the homoskedasticity 
in the first stage of the design. Figure [2] reports mean bias and rmse results for the het- 
eroskedastic design. Here cqiv-qr outperforms cqiv-ols and cqiv-dr at every quantile, which 
is expected because cqiv-ols and cqiv-dr are both misspecified for the control variable. Cqiv- 
dr has lower bias than cqiv-ols because it uses a more fiexible specification for the control 
variable. Cqiv-qr also outperforms all other quantile estimators. Most importantly, at every 
quantile, cqiv-qr outperforms both tobit estimators, which are no longer consistent given 
the heteroskedasticity in the design of the first stage. In summary, CQIV performs well 
relative to tobit in a model that satisfies the parametric assumptions required for tobit-iv to 
be efficient, and it outperforms tobit in a model with heteroskedasticy. 
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4. Empirical Application: Engel Curve Estimation 



In this section, we apply the CQIV estimator to the estimation of Engel curves. The 
Engel curve relationship describes how a household's demand for a commodity changes as the 
household's expenditure increases. Lewbel (2006) provides a recent survey of the extensive 
literature on Engel curve estimation. For comparability to the recent studies, we use data 
from the 1995 U.K. Family Expenditure Survey (FES) as in Blundell, Chen, and Kristensen 
(2007) and Imbens and Newey (2009). Following Blundell, Chen, and Kristensen (2007), 
we restrict the sample to 1,655 married or cohabitating couples with two or fewer children, 
in which the head of household is employed and between the ages of 20 and 55. The FES 
collects data on household expenditure for different categories of commodities. We focus on 
estimation of the Engel curve relationship for the alcohol category because 16% of families 
in our data report zero expenditure on alcohol. Although zero expenditure on alcohol arises 
as a corner solution outcome, and not from bottom coding, both types of censoring motivate 
the use of censored estimators such as CQIV. 

Endogeneity in the estimation of Engel curves arises because the decision to consume a 
particular category of commodity may occur simultaneously with the allocation of income 
between consumption and savings. Following the literature, we rely on a two-stage budgeting 
argument to justify the use of labor income as an instrument for expenditure. Specifically, 
we estimate a quantile regression model in the first stage, where the logarithm of total 
expenditure, is a function of the logarithm of gross earnings of the head of the household, 
Z, and demographic household characteristics, W . The control variable, V , is obtained 
using the CQIV-QR estimator in (12. lip , where the integral is approximated by a grid of 100 
quantiles. For comparison, we also obtained control variable estimates using least squares 
and probit distribution regression. We do not report these comparison estimates because 
the correlation between the different control variable estimates was virtually 1, and all the 
methods resulted in very similar estimates in the second stage. 

In the second stage we focus on the following quantile specification for Engel curve esti- 
mation: 



where Y is the observed share of total expenditure on alcohol censored at zero, 14^ is a binary 
household demographic variable that indicates whether the family has any children, and V 
is the control variable. We define our binary demographic variable following Blundell, Chen 
and Kristensen (2007)i 




Demographic variables are important shifters of Engel curves. In recent literature, "shape invariant" specifi- 
cations for demographic variable have become popular. For comparison with this literature, we also estimate 



Y.i = max(A;/3o(t/^), 0), A, = (1, A, A', ^-\Vi)), A - f/(0, 1) | A, 
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To choose the specification, we rely on recent studies in Engel curve estimation. Thus, 
following Blundell, Browning, and Crawford (2003) we impose separability between the con- 
trol variable and other regressors. Hausman, Newey, and Powell (1995) and Banks, Blundell, 
and Lewbel (1997) show that the quadratic specification in log-expenditure gives a better 
fit than the linear specification used in earlier studies. In particular, Blundell, Duncan, and 
Pendakur (1998) find that the quadratic specification gives a good approximation to the 
shape of the Engel curve for alcohol. To check the robustness of the specification to the 
linearity in the control variable, we also estimate specifications that include nonlinear terms 
in the control variable. The results are very similar to the ones reported. 

Figure [3] reports the estimated coefficients u ^ /3{u) for a variety of estimators. In addition 
to reporting results for CQIV with a quantile estimate of the control variable (cqiv), as in 
the previous numerical examples, we report estimates from the censored quantile regression 
(cqr) of Chernozhukov and Hong (2002), the quantile instrumental variables estimator with 
a quantile regression estimate of the control variable (qiv) of Lee (2007), and the quantile 
regression (qr) estimator of Koenker and Bassett (1978). We also estimate a model for 
the conditional mean with the tobit-cmle of Smith and Blundell (1986) that incorporates a 
least squares estimate of the control variable. The tobit-iv algorithm implemented in Stata 
does not converge in this application. Given the level of censoring, we focus on conditional 
quantiles above the .15 quantile. 

In the panels that depict the coefficients of expenditure and its square, the importance of 
controlling for censoring is especially apparent. Comparison between the censored quantile 
estimators (cqiv and cqr), plotted with thick light lines, and the uncensored quantile estima- 
tors (qiv and qr), plotted with thin dark lines, demonstrates that the censoring attenuates 
the uncorrected estimates toward zero at most quantiles in this application. In particular, 
censoring appears very important at the lowest quantiles. Relative to the tobit-cmle estimate 
of the conditional mean, cqiv provides a richer picture of the heterogenous effects of the vari- 
ables. Comparison of the quantile estimators that account for endogeneity (cqiv and qiv), 
plotted with solid lines, and those that do not (cqr and qr), plotted with dashed lines, shows 
that endogeneity also influences the estimates, but the pattern is more difficult to interpret. 
The estimates of the coefficient of the control variable indicate that the endogeneity problem 
is more severe in the upper half of the distribution. This is consistent with a situation where 
a strong preference to consume alcohol raises total household expenditure. 

Our quadratic quantile model is flexible in that it permits the expenditure elasticities to 
vary across quantiles of the alcohol share and across the level of total expenditure. These 



an unrestricted version of shape invariant specification in which we inchide a term for the interaction be- 
tween the logarithm of expenditure and our demographic variable. The results from the shape invariant 
specification are qualitatively similar but less precise than the ones reported in this application. 
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quantile elasticities are related to the coefficients of the model by 

ddQviu I x) = l{x'/3o(n) > 0}{/3oi(u) + 2/3o2(«) d}, 

where /3oi(w) and /3o2(m) are the coefficients of D and D^, respectively. Figure H] reports point 
and interval estimates of average quantile elasticities as a function of the quantile index u, 
i.e., u (-)■ EplddQviu \ X)]- Here we see that accounting for endogeneity and censoring 
also has important consequences for these economically relevant quantities. The difference 
between the estimates is more pronounced along the endogeneity dimension than it is along 
the censoring dimension. The right panel plots 95% pointwise confidence intervals for the 
cqiv quantile elasticity estimates obtained by the weighted bootstrap method described in 
Section 3 with standard exponential weights and B = 200 repetitions. Here we can see 
that there is significant heterogeneity in the expenditure elasticity across quantiles. Thus, 
alcohol passes from being a normal good for low quantiles to being an inferior good for high 
quantiles. This heterogeneity is missed by conventional mean estimates of the elasticity. 

In Figure [5] we report families of Engel curves based on the cqiv coefficient estimates. We 
predict the value of the alcohol share, Y, for a grid of values of log expenditure using the 
cqiv coefficients at each quartile. The subfigures depict the Engel curves for each quartile of 
the empirical values of the control variable, for individuals with and without kids, that is 

d ^ max{(l, d, rf^ w, <l>^\v)yp{u), 0} 

for {w, ), u) evaluated at w G {0, 1}, the quartiles of V for t>, and u G {0.25, 0.50, 0.75}. 
Here we can see that controlling for censoring has an important effect on the shape of the 
Engel curves even at the median. The families of Engel curves are fairly robust to the values 
of the control variable, but the effect of children on alcohol shares is more pronounced. The 
presence of children in the household produces a downward shift in the Engel curves at all 
the levels of log-expenditure considered. 



5. Conclusion 

In this paper, we develop a new censored quantile instrumental variable estimator that 
incorporates endogenous regressors using a control variable approach. Censoring and en- 
dogeneity abound in empirical work, making the new estimator a valuable addition to the 
applied econometrician's toolkit. For example, Kowalski (2009) uses this estimator to ana- 
lyze the price elasticity of expenditure on medical care across the quantiles of the expenditure 
distribution, where censoring arises because of the decision to consume zero care and en- 
dogeneity arises because marginal prices explicitly depend on expenditure. Since the new 
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estimator can be implemented using standard statistical software, it should prove useful to 
applied researchers in many applications. 

Appendix A. Notation 

In what follows ^ and 7 denote generic values for the control function and the parameter 
of the selector l(S'j'7 > It is convenient also to introduce some additional notation, 
which will be extensively used in the proofs. Let Xi{{}) := x{Di, Wi, 'd{Di^ Wi, Zi)), Si{^) := 
s{Xi{i^),Ci), := d^x{Di,Wi,v)\^=./,(^D„w„zo^ and 5*^(^9) := d^s{Xi{v),Ci)\y=^(^D,^w„Zi)- 

When the previous functions are evaluated at the true values we use Xj = Xi{{}Q), Si = 
Sii^o), Xi = Xi{^o), and Si = Si{'&o)- Also, let cpuiz) := [l{z < 0) - u]. Recall that 
A := {Y,D,W,Z,C,X,V). For a function f : A ^ R, we use ||/||oo = sup„g_4 |/(a) |; for a 
i^- vector of functions f : A ^ M.^ , we use ||/||oo = sup„g_4 ||/(a)||2. We make functions in 
T as well as estimates ^9 to take values in [0, 1], the support of the control variable V. This 
allows us to simplify notation in what follows. We also adopt the standard notation in the 
empirical process literature (see, e.g., van der Vaart, 1998), 

n 

EM]=^n[f{A)]=n-'Y,f{A,), 

i=l 

and 

n 

GM] = = n~'/' J2{fiA) - Ep[/(A)]). 

When the function / is estimated, the notation should interpreted as: 

Gn[f]=Gn[f] \f=r 

Appendix B. Proof of Theorems 1 and 2 

B.l. Proof of Theorem 1. Step 1. This step shows that ^/n0{u) - (3o{u)) = Op(l). 
By Assumption H] on the selector, for large enough n: 

l{S{dy^ > c^} < l{5'7o > - ej < l{5'7o > ^^72} < l{X'(3o{u) > C + e'} =: x, 
P-a.e., since 

\S{dy^ - S'^ol < en := Ls{\\d-Moo + \\l - 70II2) 0, 
where L5 := (||(9^s||oo V ||s||oo) is a finite constant by assumption. 
Hence, with probability approaching one, 

^(M) = arg min E„[p„(F - A(^)'/3)l(^(^9)'7 > 

flglRdim(X) 
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Due to convexity of the objective function, it suffices to show that for any e > there exists 
a finite positive constant such that 



> > 1 - e, 



where 



Let 



hm inf P inf y/nrj'En fr,,B, 

.!|r;||2 = l L 



(B.l) 



fiA) := ip^ {Y - X'/3o(n)} Xl{5'7o > 
Then uniformly in ||?7||2 = 1, 

= (1) V'^nlf] + Op(l) + r]'^/nEp[fr,^B^ 

=(2) V'^nlf] + op(l) + r]'J{u)r]B^ + T]'Qn[g\ + Op(l) 
=(3) Op(l) + op(l) + T]'J{u)r]B, + Op(l) + op(l), 

where relations (1) and (2) follow by Lemma [1] and Lemma [2] with /3 = Po{u) + B^r|/^Jn, 
respectively, using that — "i^Hoo = ov{l/^Jn)^ -i? G T, H"*? — 'i^olloo = Of>{l/^Jn) and ||/3o('w) + 
Bj]l\fn — /3o(^)||2 = 0{\l\fn)\ relation (3) holds by Chebyshev inequality. Since J{u) is 
positive definite, with minimal eigenvalue bounded away from zero, the inequality ( IB. II) 
follows by choosing B^ as a sufficiently large constant. 

Step 2. In this step we show the main result. From the subgradient characterization of 
the solution to the quantile regression problem we have 



/ 



Pnib < dim(X) max \X.-^2l^= Op(l), 

l<i<n 



fB.2) 



where 



Therefore 

Opfl) = 



/ 



/ 



f 



=(1) 



G„[/] + op(l) + v^Ep 



/ 



= (2) 



[/] + op(l) + J{u)y^{P{u) - (5o{u)) + Gnlg] + Op(l) 



where relations (1) and (2) follow by Lemma [T] and Lemma [2] with (3 = f3{u), respectively, 
using that 11^9 - = op(l/v/^), ^ e T, ||^ - = 0^(1/,/^) and - /3o{u)h = 

Op(l/v^). 
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Therefore by invertibility of J{u), 

MPi^^) - /3oH) = -J{u)-'Gn{f + 9) + op(l). 
By the Central Limit Theorem, G„(/ + g) -^d A^(0, Vaipi^f + g)), so that 
V^0{u) - Poiu)) iV(0, J{u)-WaTp{f + g)J{u)-'). 

□ 

B.2. Proof of Theorem 2. Step 1. This step shows that y/n((3^{u) — (3o{u)) = Op(l) under 
the unconditional probability P. 

By Assumption m on the selector, with probability approaching one, 

^'^(n) = arg min E„[ep„(r - X(^7/3)l(5(r)'7 > 

where e is the random variable used in the weighted bootstrap, and x = ^{^'(^oiu) > C + e'). 
Due to convexity of the objective function, it suffices to show that for any e > there exists 
a finite positive constant such that 

liminf P ( inf v^r/'E„ [ f% 1 > ) > 1 - e, (B.3) 

n^oo V||7?||2 = l L J 

where 



Let 

r{A) ■.= e-^u{Y- X'/3o(m)} Xl{5'7o > 
Then uniformly in ||?7||2 = 1, 

=(1) VG„[r] + op(l) + ti'V^^pIIIb} 

= {2) VG„[r] + 0P(1) + V'J{u)vBe + V'Gn[g'] + Op(l) 

=(3) Op(l) + op(l) + ri'J{u)riB, + Op(l) + op(l), 

where relations (1) and (2) follow by Lemma [1] and Lemma [2] with /3 = Po{u) + B^r|/^/n, 
respectively, using that — 'i9'^||oo = op(1/y^), ^'^ G T, — i^oHoo = 0]f>{l/y^) and 
||/3o(m) + B^rj/ y/n — /3o(^)||2 = 0{l/y/n)] relation (3) holds by Chebyshev inequality. Since 
J{u) is positive definite, with minimal eigenvalue bounded away from zero, the inequality 
( 1B.3[) follows by choosing B^ as a sufficiently large constant. 

Step 2. In this step we show that y^{/3'^{u) — /3q{u)) = —J{uy^Gn{f'^ + g'^) + op(l) under 
the unconditional probability P. 
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From the subgradient characterization of the solution to the quantile regression problem 
we have 



r 



Sn, W^nh < dim(X) max \\eiXi\\2/Vn = op(l) 

l<i<n 



(B.4) 



where 



Therefore 

Op(l) = ^/nEn 



r 



r 



r 



=(1) G„[r] + op(i) + v^Ep [r 

= (2) 'Gsnin + Op(l) + J{U)^C(5\U) - /3o(m)) + + Op(l), 

where relations (1) and (2) follow by Lemma [1] and Lemma [2] with (3 = I3^{u), respectively, 
using that \\d^-d^\\oo = op(l/v^), G T, \\d'-Moo = 0^(1/^) and \\^%u) - /3o{u)\\2 = 
Op(l/v^). 

Therefore by invertibility of J{u), 

V^^'iu) - /3o(m)) = - J(M)-iG„(r + g') + op(l). 

Step 3. In this final step we establish the behavior of y/n{[5^{u) — P{u)) under P^. Note 
that denotes the conditional probability measure, namely the probability measure induced 
by draws of ei, e„ conditional on the data Ai, A^. By Step 2 of the proof of Theorem 
1 and Step 2 of this proof, we have that under P: 

^0^{u)-P,{u)) = -J(M)-iG„(r+/)+op(l), ^0{u)-p,{u)) = -J{u)-^Gn{f+g)+o^{l). 
Hence, under P 

^(P^{u)-P{u)) = -J(n)-iG„(r-/+/-(?)+r„ = - J(n)-iG„((e-l)(/+(?))+r„, r„ = op(l). 
Note that it is also true that 

r„ = ope(l) in P-probability, 

where the latter statement means that for every e > 0, P'^(||r„||2 > e) — )-p 0. Indeed, this 
follows from Markov inequality and by 



Ep[r(||r„||2>e)]=P(||r„||2>e) = o(l), 
where the latter holds by the Law of Iterated Expectations and r„ = op(l). 
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By the Conditional Multiplier Central Limit Theorem, e.g., Lemma 2.9.5 in van der Vaart 
and Wellner (1996), we have that conditional on the data Ai, ...,An 

G„((e - 1) (/ + (?)) Z := iV(0, Varp(/ + ^)), in P-probability, 

where the statement means that for each z e M^^^^^) 

P'^(G„((e -!)(/ + g)) < z) Pr(Z < z). 

Conclude that conditional on the data Ai, y4„ 

Vn(i3\u) - ^{u)) ->d A^(0, J~^(M)Varp(/ + g)J-^{u)), in P-probability, 

where the statement means that for each z G R'^™*^"^) 

P"(v^(^(m) - < z) Vi{-J-\u)Z < z). 

□ 



B.3. Lemma on Stochastic Equicontinuity. 

Lemma 1 (Stochastic equicontinuity). Let e > be a positive random variable with unit 
mean and finite variance that is independent of {Y, D,W, Z, X,V) , including as a special 
case e = l, and set, for A = (e, Y, D, W, Z, X, V) and x = 1(^'/3o(m) > C + e'), 

/(A ^9, /3, 7) := e ■ [1{Y < X(^)'/3) - u] ■ X(^) • 1(S(^)'7 > ■ X- 

Under the assumptions of the paper, the following relations are true. 

(a) Consider the set of functions 

T = {f{A, ^, /3, 7)'« : (^9, /3) G To X i3, 7 e r, « e M"'"^^^), H^lh < 1}, 

where T is an open neighborhood of '-/q under the \\ ■ II2 metric, B is an open neigh- 
borhood of Pq{u) under the || ■ II2 metric, Tq is the intersection of T, defined in 
Assumptions^ with a small neighborhood of under the || ■ ||oo metric, which are 
chosen to be small enough so that: 

\X{'d)'P - X'Pq{u)\ < e/2, P-a.e. V(?9, /3) G Tq x i3, 

where e' is defined in Assumption 4- This class is P-Donsker with a square integrable 
envelope of the form e times a constant. 

(b) Moreover, if (t?, /3,7) — )■ (-i^o, /3o(w), 70) in the \\ ■ ||oo V || ■ II2 V || ■ II2 metric, then 

WfiA, ^, /3, 7) - fiA, ^?o, /9o(w), 7o) \\p,2 ^ 0. 
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(c) Hence for any (3, 7) — j-p {{}q, /3o(m), 70) in the || • ||oo V || ■ II2 V || ■ II2 metric such that 
^ G To , 

||G„/(A ^, (3, 7) - Gnf{A, ^0, (3o{u),7o) h 0. 

(d) For for any {^,(3,^) — )-p {-do, (3o{u),'yo) in the \\ ■ ||oo V || ■ II2 V || ■ II2 metric, so that 



\W — ^'^Woo = op{l/^/n), where -de To, 

we have that 

||G„/(A^,A7) -Gn/(A^o,/3o(w),7o)||2 ^pO. 



Proof of Lemma [T]. The proof is divided in proofs of each of the claims. 
Proof of Claim (a). The proof proceeds in several steps. 
Step 1. Here we bound the bracketing entropy for 

Xi = {[1{Y < X(^)'/3) -u]x:^eB,^e To}. 

For this purpose consider a mesh {"^k} over To of || ■ ||oo width 6, and a mesh {/3i} over B of 
II • II 2 width 6. A generic bracket over Xi takes the form 

[^?, z\] = [{1{Y < - ^6) - u}x, {1{Y < + ^S) - u}x], 

where k = Lx max^gg ||/3||2 + Lx, and Lx '■= ||5^a;||oo V ||a;||oo- 

Note that this is a valid bracket for all elements of Xi induced by any located within 5 
from i)k and any /3 located within 6 from Pi, since 

|X(^)'/3-X(i9fc)'A| < \{X{^)-X{^k)yf3\ + \X{^ky{/3-/3k)\ 

< LxSmax\\m2 + Lx6 < k6, (B.5) 

and the L2{P) size of this bracket is given by 

\K-tl\\p,2 < VMP{y e [Xi^kYPi ± k6] I D, W, Z,C,x = 1}] 

< JM sup P{Y e[y±^6]\X,C,x=l}] 

Y s/e(C+K5,oo) 

< V||/y(- I •)l|oo2/€5, 



provided that 2k5 < e'/2. In order to derive this bound we use the condition |X(79)'/3 — 
X'f3o{u)\ < e'/2, P-a.e. W{^,(3) E Tq x B, so that conditional on x = 1 we have that 
X{^)'/3 >C + e'/2; and 

P{Y e-\D,W,Z,C,x=l} = P{Y e-\D,w,z,v,C,x = i} = P{Y e-\x,C,x = i}, 
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because V — ■&q{D, W, Z) and the exclusion restriction for Z. Hence, conditional on X, C and 
X = 1, y does not have point mass in the region [X{'dk)' j3i±K8\ C (C, oo), and by assumption 
the density of Y conditional on X, C is uniformly bounded over the region (C, oo). 

It follows that 

logiVo(e,Xi,L2(P)) < logiV(e2,To, || ■ |U) + log iV(e, .B, || ■ II2) < l/le^ log^ e) + log(l/e), 

and so Xi is P-Donsker with a constant envelope. 
Step 2. Similarly to Step 1, it follows that 

X2 = {X{:&)'a : ^? e To, a e M'^^'"^, ||a||2 < 1} 

also obeys a similar bracketing entropy bound 

log7Vo(^'^2,i^2(P)) < 1/(6^^6) +log(l/e) 

with a generic bracket taking the form [«2?'^2] — [^{'^k)' Pi — 1^^, ^i'^k)' Pi + i^^]- Hence, this 
class is also P-Donsker with a constant envelope. 

Step 3. Here we bound the bracketing entropy for 

X3 = {iWy7>^) :^eTo,7er}. 

For this purpose consider the mesh {'&k} over Tq of || ■ ||oo width 6, and a mesh {7;} over F 
of II ■ II 2 width 5. A generic bracket over X3 takes the form 

[I'^il] = [i(5(^fe)'7i - > 0, HSii^kYii + Kd> 0], 

where k — Ls max^gr II7II2 + Lg, and L5 := ||9^s||oo V ||s||oo- 

Note that this is a vahd bracket for all elements of X3 induced by any ■& located within 5 
from and any 7 located within 5 from 7^, since 

\s{d)'^-s{§k)'ii\ < \{s{^)-si'dk)yi\ + \sW{7-7k)\ 

< L55 max II7II 2 + < AC(^, (B.6) 
and the size of this bracket is given by 

K - il\\p,2 < VP{\SWli-<i\<2KS} < ^fs2KS, 

where fs is a constant representing the uniform upper bound on the density of random 
variable 5'('j?)'7, where the uniformity is over 1? e Tq and 7 e F. 

It follows that 

logN[]{e,Is,L2{P)) < logA^(e^To, || • lU) + log Ar(e, F, || • II2) < l/(e2 V e) + log(l/e) 
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and so X3 is P-Donsker with a constant envelope. 

Step 4. In this step we verify the claim (a). Note that J-" = e ■ Xi ■X2 -X^. This class has a 
square-integrable envelope under P. 

If we make a further assumption that the random variable e obeys Ep|ep"'"'^ < 00, for 
5 > 0, then the class J-" is P-Donsker by the following argument. Note that the product 
Xi ■ X2 ■ X3 of uniformly bounded classes is P-Donsker, e.g., by Theorem 2.10.6 of van der 
Vaart and Wellner (1996). Under the stated assumption the final product of the random 
variable e with the P-Donsker class remains to be P-Donsker by the Multiplier Donsker 
Theorem, namely Theorem 2.9.2 in van der Vaart and Wellner (1996). 

It is also interesting to give a direct argument based on bracketing, which also requires the 
somewhat weaker assumption of e > having finite variance. Indeed, our generic brackets 
over J-" take the form 



mm e ii, max e in ii 
(9Ai)e{o,i}3 (^,/^,i)e{o,l}3 



For example, if i5;'^25^3 — O5 then = e i\ i\ i\ and = e i\ i\ i\. Note that these are 
valid brackets for all / that are induced by an element of Xi located in the bracket [^5,^1], 
an element of X2 located in the bracket [i^?^]' ^^"^ element of X3 located in the bracket 
[ig, ig]. Using that the constructed brackets are uniformly bounded, by say A;, the size of this 
bracket is 

- /1p,2 < V^kW\i\ - zj||p,2 + ¥l - i\\\p,2 + ¥l - 4||P,2) < v^E^^fc^s^, 

where we used the independence of random variable e from random vector A. Therefore, 
logiV[](e,^,L2(P)) <logiV[](e,Xi,L2(P)) + logiVo(e,X2,L2(P)) + logiV[](e,X3,L2(P)), 

and the class is thus P-Donsker by the previous steps. 

Proof of Claim (b). The claim follows by the Dominated Convergence Theorem, since any 
/ G J-" is dominated by a square-integrable envelope under P, and by the following three 
facts: 



(1) in view of the relation such as (iRSj) . 1{Y < X{^yf3)x -> 1{Y < X'f3o{u))x every- 



where, except for the set {A ^ A : Y = X'/3o(m)} whose measure under P is zero by 
Y having a uniformly bounded density conditional on X,C; 

(2) in view of the relation such as (IB.Sp . |X('i9)'/3 — X'l3o{u) \ — )■ everywhere; 

(3) in view of the relation such as (]B.6p . l(S'('t9)'7 > <j) — )■ l(5"7o > ?) everywhere, 
except for the set {A ^ A : S"7o = <j} whose measure under P is zero by S"7o having 
a bounded density. 
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Proof of Claim (c) . This claim follows from the asymptotic equicontinuity of the empirical 
process (G„[/], f & T) under the I/2(-P) metric, and hence also with respect to the || • ||oo V 
II • II2 V II • II2 metric in view of Claim (b). 

Proof of Claim (d). It is convenient to set / := /(A, /3, 7) and / := /(^, /3, 7). Note 
that 

|Gn[/-/]| < |V^E„[/-7]| + |y^Ep(/-7)| (B.7) 

< v^E„[C] + y^Ep[C] (B.8) 

< G„[C]+2V^Ep[C], (B.9) 

where |/| denote an application of absolute value to each element of the vector /, and ( is 
defined by the following relationship, which holds with with probability approaching one, 

\f-f\<\e\-\\X{d)-X{d)y + g + h<C:=e-LxAn + g + h, A„ > ||^- ^|U, (B.IO) 

where Lx — ||9^a^||oo V ||a;||oo) and, for some constant k, 

g-e- 1{\Y - X{^yp\ < fcAJx, and h -.^ e ■ l{\S{dy^ - ^| < A;A„}, 

and A„ = o{l/y/n) is a deterministic sequence. 

Hence it suffices to show that the result follows from 

Gn[C ] = op(l), (B.ll) 

and 

V^Ep[C] =op(l). (B.12) 

Note that since A„ 0, with probability approaching one, ^ and h are elements of the 
function classes 

g = {e.l(|y-X(^y/3|<A;)x:^eTo,/3eB,A;e[0,e74]}, 

n = {e-ii\S{^)'j-q\<k):'&eTo,jer,ke[0,i]}. 

By the argument similar to that in the proof of claim (a), we have that 

log7V[](e,^,L2(P)) < l/{eHog'e) and log NQ{e,H, L^iP)) < l/{eHog'e). 
Hence these classes are P-Donsker with unit envelopes. Also note that 



||?||p,2 < VEp[e2] • P{\Y - < kAn} < V^LMT^U^ = o(l), (B.13) 

||^||p,2 < -v/Ep[e2] • P{\S0yj-<^\ < kAn} < ^^fskAn = o(l), (B.14) 
by the assumption on bounded densities and Ep[e^] = 2. 
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Conclude that the relation fiRTTD holds by flRTOjl . (lRl3|) . fiRlil) . the P-Donskerity of the 
empirical processes ((G-„[/i], h El-i) and (G„[(?], (7 G ^) and hence their asymptotic equiconti- 
nuity under the || ■ \\p^2 metric. Indeed, these properties imply 



|P,2 = 0(1) ^ Gn[C ] = 0(1). 

To show (]B.12p note that 

\/^Ep|C| < Ep|e| ■ Lx\fn^n + ||?||p,i + 
since A„ = o{l/y/n), and 

||?||p,i < Ep|e| ■P{\Y- X{d)'l3\ < fcA„} < 2fc||/yi 
||^||p,i < Ep|e| ■ P{\S(d)'^ - ^1 < fcAJ < 2kfs/^n 
by the assumption on bounded densities. 



P,i 



0(1), 



1 00 An 



:i/v^), 



o(l/v^) 



B.4. Lemma on Local Expansion. 

Lemma 2 (Local expansion). Under the assumptions stated in the paper, for 
5 = - /3o{u)) =Op(l);7 = 7o + op(l); 

A{d,w,z) = y/^{d{d,w,z)-^Q{d,w,z)) = ^ E„[£(A f^, z)] + op(l) m r{Vll), 
l|V^E„[£(A,-)]||oo = Op(l), 

vje have that 

^ Epl{r < X0)'P}X0)l{S0)'-i > <;}x = J{u)5 + ^ E„ [g{A)] + op(l), 

where 

g{A)= j B{a)e{A,d,r)dP{a,d,r), 5(A) := /y(X%(M)|X, C)XX%(m)1(^'7o > 0- 

Proof of Lemma [2]. We have that with probability approaching one, 

1{^(^) 7 > ^ < HS'lo > ^ - ej < l{5'7o > <i/2} < l{X'/3o(w) > C + e'} = x, 
P-a.e., by assumption on the selector since 

|5(^)'7 - 5'7o| < en := Ls{0 - Moo + II7 - 70II2) 0, 
where L5 = (||9t,s||oo V ||s||oo) is a finite constant by assumption. 
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Hence uniformly in X such that X' I3q{u) > C + e', 
x/nEp[l{r < X{d)'P} I D, W, Z, C] 

= fY{X{dxyPx I D,W,Z,C){X{dx)'5 + X{dx)'PxHD,W,Z)} 
= fY{X'/3o{u) I D, W, Z, C){X'6 + rPo{u)A{D, W,Z)} + Rx, 
= /y (X'/3o(m) I X, C){X'6 + X'/3o{u)A{D, W, Z)} + Rx, 
R= sup \Rx\=op{l), 

{X:X'l3o{u)>C+€'} 

where -t^x is on the hne connecting and i) and Px is on the hne connecting Po{u) and /3. 
The first equahty follows by the mean value expansion. The second equality follows by the 
uniform continuity assumption of /y(- | X, C) uniformly in X, C, uniform continuity of X{-) 
and X{-), and by — i^oHoo — >p and \\/3 — /3q{u)\\2 — J-p 0. The third equality follows by 

/y(- I D,W,Z,C) = fY{- I D,W,Z,V,C) = fy{-\X,C) 

because V = 'do{D, W, Z) and the exclusion restriction for Z. 

Since /y (■ | ■) and the entries of X and X are bounded, and 6 = Op{l) and || A||oo = Op(l), 
with probability approaching one 

Epi(r < x{^yp)x{^)i{s{^y^ > 

= Ep[fY{X'f3o{u)\X,C)XX'l{S{dyj > c;}]6 

+Ep[fY{X'Po{u)\X,C)XX'Po{u)l{S{dy^ > <;}A{D,W,Z)]+Of{R). (B.15) 

Furthermore since 

Ep|l(S'7o > ^) - l{S{dy^ > 01 < Ep[l(5'7o e ± e„])] < fstn 0, 

where fs is a constant representing the uniform upper bound on the density of random 
variable S"7o, the expression flB.lSp is equal to 

J{u)6 + Ep[/y(X'/3o(M)|X, C)XX'Po{u)l{S'^o > OA(D, W, Z)] + Op{fsen + R). 

Substituting in A{d,w,z) = -Jn E„[£(y4, d, w, 2;)] + Op(l) and interchanging Ep and E„, we 
obtain 

Ep[/y(X'/3o(M)|X,C)XX'/3oHl(S'7o > OA(D,W^,Z)] = E„[(?(A)] + op(l). 
The claim of the lemma follows. □ 



Appendix C. Proof of Theorem 3 



35 



To show claim (1), we first note that by Chernozhukov, Fernandez- Val and Melly (2009), 

V^{n{v) - TToiv)) = ^ Ve.Ep [foiR'M''^) I R)RRT' b - < R'M'")}]R- 

By the Hadarmard differentiabihty of rearrangement-related operators in Chernozhukov, 
Fernandez- Val and Gahchon (2010), the mapping tt ^ 0^ from £'^{0,1)'^^^^^ to l°°{Vn) 
defined by 

0^(d,r) = / l{rV(T;) < d}dv 

J(o,i) 

is Hadamard difTcrcntiable at tt = ttq, tangentially to the set of continuous directions, with 
the derivative given by 

<P^o[h] = -fD{d\ryh{Md,r)), 

where '&Q{d,r) = 0^g((i, r). Therefore by the Functional Delta Method (Theorem 20.8 in van 
der Vaart, 1998), we have that in £°°(D7^), for d{d,r) = (f)^{d,r), 

V^{d{d, r) - Md, r)) = -fD{d \ r)r'^ EILi e.Ep [foiR'noiMd, r)) | R)RRr' X 

x[do{d,r) - 1{D < R'no{Md,r))}]R + op{l). 

The claim (1) then follows immediately. Also for future reference, note that the result also 
implies that 

x/^(7r(-) -7ro(-)) ^ ^vr in £°°(0,1), and r'V^(7r(-) - 7ro(-)) ^ rX in £°°(7^ x (0, 1)), 

(C.l) 

where is a Gaussian process with continuous sample paths. 
The proof of claim (2) is divided in several steps: 

Step 1. In this step we construct T and bound its entropy. Let C|^(0, 1) denote the class 
of functions / : (0, 1) — ?> M with all derivatives up to order 2 bounded by a constant M, 
including the zero order derivative. The covering entropy of this class is known to obey 
logiV(e, Clio, 1), II ■ Hoc) < e-i/^ Hence also logiV(e, xf^^^'^CUO, 1), || • ||oo) < e"'/'- Next 
construct for some small m > 0: 

^ " {/o 1) ^^^'''^''^ - ^^"^"^ ■ " "'"'^^''^^ ^ X^t"/'')C^(0, 1), R'dn{v) > m P-a.e.| 

Then for any tt and tt obeying the conditions in the display and such that ||7r — 7f ||oo < 5, 
/ 1{R!'k{v) < D}dv - I 1{R't^{v) < D}dv 

i(0,l) -'(0,1) 

< / llR'niv) -De [-WRhS, \\R\\2S]}dv < -{{RhS < 5, 

J(Q,1} ^ J ^ 



36 

P-a.e., since the density of r'7r{V), V ~ C/(0, 1), is bounded above by 1/m. We conclude 
that 

log7V(e, T, II . lU) < log7V(e, xjr/^)cl,(0, 1), || • |U) < e'^^. 

Step 2. In this step we show that there exists E T such that ||'J? — 'J?||oo = op(l/\/n)- 
We first construct n such that 

^A^||7f - 7r||oo = op(l), and max v^||r'(7f - 7f)||oo = op(l), (C.2) 

where with probabihty approaching one, % e x5r/^)C|^(0, 1) and R'd%{v) > m P-a.e., for 
some M and some m > 0. 

Wc construct vf by smoothing tt component by component. Let the components of vr 
be indexed by 1 < j < dim(i?). Before smoothing, we need to extend tt^ outside (0, 1). 
Start by extending the estimand ttoj outside (0, 1) onto the e-expansion (0, 1)^ smoothly 
so that the extended function is in the class C^. This is possible by first extending d'^Tioj 
smoothly and then integrating up to obtain lower order derivatives and the function. Then we 
extend the estimator ttj to the outer region by simply setting Tfj{v) = 7Toj{v) over v ^ (0, 1). 
Note that the extended empirical process y/n{7rj{v) — noj{v)) remains to be stochastically 
equicontinuous, since it vanishes asymptotically at the boundary of (0, 1). Then we construct 
TTj as the smoothed version of tTj, namely 

TTjiv) = / Tfj{z)[K{{z - v)/h)/h]dz, v e (0, 1), 

where < h < e is bandwidth such that ^Jnlv' — )■ and ^Jn}^ — t-oo; i^:M— )-Misa third 
order kernel with the properties: d^K are continuous on [—1, 1] and vanish outside of [—1, 1] 
for ^ = 0, 1, 2, J K{z)dz = 1, and J z^K{z)dz = for = 1, 2. Such kernel exists and can be 
obtained by reproducing kernel Hilbert space methods or via twicing kernel transformations 
(Berlinet, 1993, and Newey, Hsieh, and Robins, 2004). We then have 

Vn{7rj{v) - TTjiv)) = / Vn[7^j{z) - 7roj{z) - {nj{v) - 7rojiv))][K{[z - v]/h)/h]dz 

+ / Vn{noj{z) - noj{v))[K{[z - v]/h)/h]dz. 

The first term is bounded uniformly in v G (0, 1) by a;(2/i)||i^||oo where 
uj{2h) = sup \^/n[^fj{z) - 7roj{z) - {7fj{v) - 7roj{v))] \ ^p 0. 

\z-u\<2h 

The second term is bounded uniformly in v e (0, 1), up to a constant, by 

Vn\\d^noj\\ooh^ I X^K{X)d\ < ^/^h? 0. 
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This establishes the equivalence (lC.2p . in view of compactness of TZ. 

Next we show that ||c?^7fj||oo < 2||(9^7roj||oo =: M with probability approaching 1. Note 
that 

d%{v) - d^noj{v) = [ nj{z)[d''K{[z - v]/h)/h^]dz - d'^nojiv), 
which can be decomposed into two pieces: 

n-'^'h~' /(o,i). nV2(vf,.(^) _ rro,{z))[d'K{[z - v]/h)/h]dz 
+ J^^^^y[d'K{[z - v]/h)/h^]no,iz)dz - d'no.iv). 

The first piece is bounded uniformly in f G (0,1) by n~^/^h~^uj{2h)\\d'^K\\^ -> 0, while, 
using the integration by parts, the second piece is equal to 

/ [d^nojiz] - d\oA^)][Ki[z - v]/h)/h]dz. 
This expression is bounded in absolute value by 

\\K\\^ sup \d'^TToj{z) - d'^TTojiv)\ -> 0, 

\z~v\<2h 

by continuity of d'^nQj and compactness of (0, l)^ Thus, we conclude that ||<9^7fj— 9^vroj||oo — >p 

0, and we can also deduce similarly that Wdnj — SttojUoo — 0, all uniformly in 1 < j < 
dim(i?), since dim(i?) is finite and fixed. 

Finally, since by Assumption lU^b) the conditional density is uniformly bounded above by 
a constant, this implies that R'dirolv) > k P-a.e., for some constant k > 0, and therefore we 
also have that with probability approaching one, R'dTf{v) > k/2 =: m P-a.e. 

Next we construct 

^(d,r) = (j)^{d,r) = / l{r'TT{v) < d}dv. 

J {0,1) 

Note that by construction -t? G T for some M with probability approaching one. It remains 
to show the first order equivalence with -(9. 

By the Hadarmard differentiability for the mapping (p^^ stated earlier and by the functional 
delta method (Theorem 20.8 in van der Vaart, 1998), and have the same first order 
representation in i'^{V7l), 

v^m - M-)) = Mk-) - M-)) + op(i), 

1. e., v^||^9-?||oo ^pO. □ 
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Appendix D. Proof of Theorem 4 

Claim (1) follows from the results of Chernozhukov, Fernandez- Val and Melly (2009). Also 
for future reference, note that these results also imply that 

VE{7r{-)-7ro{-))^ Z^m£°°{V), and r'V^(7r(-) - 7ro(-)) ^ r% in £~(7eP), (D.l) 

where is a Gaussian process with continuous sample paths. 
The proof of claim (2) is divided in several steps: 

Step 1. In this step we construct T and bound its entropy. Let C£f(P) denote the class 
of functions / : — )■ M with and all the derivatives up to order 2 bounded by a constant 
M, including the zero order derivative. The covering entropy of this class is known to obey 
log(e,C2,(D),||.|U) <e-V2. Hence 

log(6, xjLr^Cl,(P),||-|U)<6-V^ 

Next construct 

T = {a{R'7t{D)) : TT = (TTi, ...,7rdi„.(«)) e xfrl^'^CliV)] . 

Then, for any tt and tt obeying the condition in the definition of the preceding class such 
that IItt — 7f lloo < S, 

\A(R'7r{D)) - A{R'ti{D))\ < \\dA\\^ sup \\r\\^S. 

reiz 

We conclude that 

log7V(e,T, II • lU) < log7V(e, x^t^^^)c|,(P), || • |U) < e'V^. 

Step 2. In this step we show that there exists e T such that Hi? — 'i?||oo — op(l/-\/ri)- 
We first construct tt and tt such that, 

y^||7f - 7r||oo = op(l), and max yn||r'(7r - 7r)||oo = op(l), (D.2) 

where with probabihty approaching one, tt G 'xf=i^^Cl[(V), for some M. 

We construct tt by smoothing tt component by component. Before smoothing, we extend 
the estimand ttq^ outside V, onto the e-expansion smoothly so that the extended function 
is of class C^. This is possible by first extending the third derivative of ttoj smoothly and 
then integrating up to obtain lower order derivatives and the function. Then we extend ifj 
to the outer region by simply setting nj{d) — 7roj{d) over d ^ T>. Note that the extended 
process y/n{Tfj{d) — Troj{d)) remains to be stochastically equicontinuous, since it vanishes 
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asymptotically at the boundary of V. Then we define the smoothed version of Tij as 

^■{d) = [ nj{z)[K{{z - d)/h)/h]dz, deV, 

where < h < e is bandwidth such that ^Jnh? — )■ and y/nh"^ — >oo; J-Misa third 

order kernel with the properties: d^K are continuous on [—1, 1] and vanish outside of [—1, 1] 
for /i = 0,1,2, J K{z)dz = 1, and J z^K{z)dz = for /i = 1,2. Such kernel exists and 
can be obtained by reproducing kernel Hilbert space methods or via twicing kernel methods 
(Berlinet, 1993, and Newey, Hsieh, and Robins, 2004). We then have 



y/n{'Kj{d) - TTjld)) = / ^/n[7^j{z) - 7Toj{z) - {7rj{d) - TToj {d))][K{[z - d]/h)/h]dz 

+ / v^(vroj(^) - TiQj{d))[K{\z - d]/h)/h]dz. 
The first term is bounded uniformly in (i G P by w(2/i)||i^||oo — ^-p where 



u{2h) = sup \^/n[nj{z) - 7roj(z) - {nj{d) - vroj((i))]| -)-p 0. 

\z-u\<2h 

The second term is bounded uniformly in d G P, up to a constant, by 

V^||5%||oo/i'y" X^K{X)dX < y/Eh' ^ 0. 

This establishes the equivalence (ID. 20 . in view of compactness of TZ. 

Next we show that ||9^7rj||oo < 2||(9^7roj||oo := M with probability approaching 1. Note 
that 



d%{d) - d'noj{d) = / 9j{z)[d'K{[z - d]/h)/h']dz - d'nojid), 
which can be decomposed into two pieces: 

n-V2/,-2 j^^ n^l^{nj{z) - 'K^j{z))[d^ K{[z - d]/h)/h]dz 
+ ^^^[d''K{\z - d]/h)/h^]'KQ,{z)dz - d^TTo^id). 

The first piece is bounded uniformly in ci G by n^^^'^h^'^uj{2h)\\d'^K\\oo — > 0, while, using 
the integration by parts, the second piece is equal to 

[d^7ioj{z) - d\ojid)][K{[z - d]/h)/h]dz, 

which converges to zero uniformly in c/ G P by the uniform continuity of d'^iToj on and 
by boundedness of the kernel function. Thus Wd'^rcj — 9^vroj||oo "~^p 0, and similarly conclude 
that \\dTTj — dnojWoo — >p 0, where convergence is uniform in 1 < j < dim(i?), since dim(_R) is 
finite and fixed. 
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We then construct d{r,d) = A{r'TT{d)). Note that by the preceding arguments -(9 G T 
for some M with probabihty approaching one. Finally, the first order equivalence \/n||'(9 — 
''^lloo — follows immediately from f lD.2|) . boundedness of ||9A||oo and compactness of TZ. 
□ 



Appendix E. Monte Carlo designs 

For the homoskedastic design, we use the following parametric linear version of the system 
of equations f l2.4p - fl2.5p to generate the observations: 

D = TToo + TToiZ + 7ro2W + <!>-\V), V^-f/(0,l), (E.l) 
Y* = f3oo + f3oiD + f3o2W + <^-\e), e^U{0,l), (E.2) 

where denotes the quantile function of the standard normal distribution, and 

^^^(e)) is jointly normal with correlation po- Though we can observe Y* in the simulated 

data, we artificially censor the data to observe 

Y = max(y*, C) = max(/3oo + f^oiD + f3o2W + <i>-\e), C). (E.3) 

From the properties of the multivariate normal distribution, $^^(e) = Pq^~^{V) + (1 — 
p2)i/2$-i(f/)^ where U - f/(0, 1). Using this expression, we can combine flE.ip and ( IE.3P 
for an alternative formulation of the censored model in which the control term Vi is included 
in the equation for the observed response: 

Y = max(r*, C) = max(/3oo + PoiD + P02W + po^-\V) + (1 - piy/^^-\U), C). 

This formulation is useful because it indicates that when we include the control term in the 
quantile function, its true coefficient is po- 

In our simulated data, we create extreme endogeneity by setting po = .9. We set ttoo = 
/3oo = 0, and ttqi = 71q2 = /3oi = /3o2 = 1- We draw the disturbances [^""^(V), ^""^(e)] from a 
bivariate normal distribution with zero means, unit variances and correlation po- We draw 
Z from a standard normal distribution, and we generate to be a log-normal random 
variable that is censored from the right at its 95th percentile, qw Formally, we draw W 
from a standard normal distribution. We then calculate qw = Qvk(-95), which differs across 
replication samples. Next, we set W = m.m{e^ , qw) ■ For comparative purposes, we set the 
amount of censoring in the dependent variable to be comparable to that in Kowalski (2009). 
Specifically, we set C = C = Qy. (.38) in each replication sample. We report results from 
1,000 simulations with n = 1, 000. 
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For the heteroskedastic design, we replace the first stage equation for D in ( ]E.1|) by the 
following equation: 

= 7roo + 7roiZ + 7ro2W^ + (7ro3 + 7ro4Vr)$"'(V), ^ - f/(0, 1) (E.4) 
where we set ttos = 7ro4 = 1. The corresponding conditional quantile function is 

Qoiy \ W,Z)= TToo + TToiZ + TToaW^ + (7ro3 + T^QiW)^^^{v), 

and can be consistently estimated by quantile regression or other estimator for location-scale 
shift models. 

Appendix F. CQIV Technical Details and Robustness Diagnostic Test 

For the OLS estimator of the control variable, we run an OLS first stage and retain 
the predicted residuals from the OLS first stage as the control variable. For the quantile 
estimator of the control variable, we run first stage quantile regressions at each quantile from 
.01 to .99 in increments of .01. Next, for each observation, we compute the fraction of the 
quantile estimates for which the predicted value of the endogenous variable is less than or 
equal to the true value of the endogenous variable. We then evaluate the standard normal 
quantile function at this value and retain the result as the estimate of the control variable. 
In this way, the quantile estimate of the control variable allows for heteroskedasticity in the 
first stage. 

For the distribution regression estimator of the control variable, we first create a matrix 
n * n oi indicators, where n is the sample size. For each value of the endogenous variable 
in the data set yj in columns, each row i gives if the log-expenditure of the individual i is 
less or equal than yj {l{yi < yj))- Second, for each column j of the matrix of indicators, 
we run a probit regression of the column on the exogenous variables. Finally, the estimate 
of the control variable for the observation i is the quantile function of the standard normal 
evaluated at the predicted value for the probability of the observation i = j. 

In Table \B1\ we present the CQIV robustness diagnostic tests suggested in section [3] for 
the CQIV estimator with an OLS estimate of the control variable. In our estimates, we 
used a probit model in the first step, and we set go = 10 and gi = 3. In practice, we do not 
necessarily recommend reporting the diagnostics in Table [BT| but we have included them here 
for expositional purposes. In the top section of the table, we present diagnostics computed 
after CQIV Step 1. At the 0.05 quantile, observations are retained in Jq if their predicted 
probability of being uncensored exceeds 1 — u + ko = 1 — . 05 + .0445 = .9945. Empirically, this 
leaves 47.0% of the total sample in Jq in the median replication sample. In all statistics, the 
variation across replication samples appears small. However, as intended by the algorithm, 
there is meaningful variation across the estimated quantiles. As the estimated quantile 
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increases, the percentage of observations retained in Jq increases. From these diagnostics, 
the CQIV estimator appears well-behaved in the sense that the percentage of observations 
retained in Jq is never very close to or 100. 

In the second section of Table IBll we present robustness test diagnostics computed after 
CQIV Step 2. Observations are retained in Ji if the predicted Yi exceeds Ci + iji, where 
the median value of Ci, as shown in the table, is 1.60, and the median value of <ji at the 
.05 quantile is 1.70. As desired, at each quantile, the percentage of observations retained in 
Ji is smaller than the percentage of observations with predicted values above Ci but larger 
than the percentage of observations retained in Jq. As shown in sections of the table labeled 
"Percent Jq in Ji" and "Count Ji not in Jq" Jq is almost a proper subset of Ji. 

In the last section of Table iBlt we report the value of the Powell objective function obtained 
after CQIV Step 2 and CQIV Step 3. The last column shows that on average the final CQIV 
step represents an improvement in the objective function in 36-51% of replication samples 
across the estimated quantiles. In our CQIV simulation results, we report the results from 
the third step. Researchers might prefer to select select results from the second or third step 
based on the value of the objective function. 
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Figure 1. Homoskedatic design: Mean bias and RMSE of Tobit and QR 
estimators. Results obtained from 1,000 samples of size n = 1, 000. 
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Figure 2. Heteroskedastic design: Mean bias and RMSE of Tobit and QR 
estimators. Results obtained from 1,000 samples of size n = 1, 000. 
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Figure 3. Coefficients of Engel Curves 
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Figure 5. Family of Engel curves: each panel plots Engel curves for the three 
quantiles of alcohol share. 
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Table B1. CQIV Robustness Diagnostic Test Results 

for CQIV with OLS Estimate of the Control Variable - Homoskedastic Design 



CQIV-OLS Step 1 







1, 

Ko 






Percent JO 










Quantile 


IVIedian 


Min 


Max 


Median 


Min 


Max 








0.05 


0.04 


0.04 


0.05 


47.20 


43.30 


50.30 








0.1 


0.09 


0.06 


0.10 


49.10 


46.00 


51.30 








0.25 


0.20 


0.15 


0.24 


52.20 


50.50 


53.70 








0.5 


0.36 


0.26 


0.46 


55.80 


54.80 


56.80 








0.75 


0.43 


0.29 


0.58 


59.40 


57.70 


61.10 








0.9 


0.37 


0.22 


0.58 


62.40 


60.30 


65.10 








0.95 


0.30 


0.18 


0.54 


64.20 


61.40 


67.50 








CQIV-OLS Step 2 




























Percent J1 




Percent Predicted Above 


Quantile 


Median 


Min 


Max 


Median 


Min 


Max 


Median 


Min 


Max 


0.05 


1.70 


1.45 


2.01 


50.70 


46.70 


54.90 


52.30 


48.20 


56.70 


0.1 


1.71 


1.44 


1.96 


52.80 


49.50 


55.50 


54.50 


51.10 


57.30 


0.25 


1.71 


1 .4b 


1 .yo 


56.30 


53.60 


oo./U 


58.10 


55.30 


60.60 


0.5 


1.72 


1.44 


2.02 


60.10 


57.60 


63.40 


62.00 


59.40 


65.40 


0.75 


1.73 


1.47 


1.99 


64.00 


61.20 


66.80 


66.00 


63.10 


68.90 


0.9 


1.75 


1.44 


2.01 


67.40 


64.60 


70.60 


69.50 


66.60 


72.80 


0.95 


1.76 


1.49 


2.02 


69.30 


65.60 


72.80 


71.50 


67.70 


75.10 






C 




Percent JO in J 1 




Count in J1 not in JO 


Quantile 


Median 


Min 


Max 


Median 


Min 


Max 


Median 


Min 


Max 


0.05 


1.60 


1.33 


1.85 


100 


97.7 


100 


36 





81 


0.1 


1.60 


1.33 


1.85 


100 


99.0 


100 


37 


7 


74 


0.25 


1.60 


1.33 


1.85 


100 


99.6 


100 


40 


15 


68 


0.5 


1.60 


1.33 


1.85 


100 


99.6 


100 


43 


23 


78 


0.75 


1.60 


1.33 


1.85 


100 


99.7 


100 


47 


17 


74 


0.9 


1.60 


1.33 


1.85 


100 


99.7 


100 


50 


15 


88 


0.95 


1.60 


1.33 


1.85 


100 


99.1 


100 


51 


16 


97 


Comparison of Objective Functions 


















Objective Step 3 




Objective Step 2 




Objective Step 3<0bjective Step 2 


Quantile 


Median 


Min 


Max 


Median 


Min 


Max 


Median 




Mean 


0.05 


5058 


4458 


5674 


5054 


4400 


5753 







0.44 


0.1 


8939 


7925 


9946 


8927 


7888 


10049 







0.47 


0.25 


17292 


15100 


19839 


17271 


14741 


20052 







0.44 


0.5 


22859 


18692 


27022 


22837 


18306 


27091 







0.45 


0.75 


16073 


9603 


22872 


15895 


8737 


22866 







0.42 


0.9 


-1016 


-9624 


7150 


-1047 


-10834 


9265 







0.45 


0.95 


-13815 


-24602 


-2884 


-14034 


-27816 


-1919 







0.44 



N=1,000, Replications=1,000 



